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We show that the neutrino models, as suggested by Low which have an additional Abelian 
family symmetry and a real Higgs singlet to the default see-saw do not hinder the possibility of 
successful thermal leptogenesis. For these models (neglecting radiative effects), we have investigated 
the situation of strong washout in both the one-flavor approximation and when flavor effects are 
included. The result is that while such models predict that #13 = and that one light neutrino to be 
massless, they do not modify or provide significant constraints on the typical leptogenesis scenario 
where the final asymmetry is dominated by the decays of the lightest right-handed neutrinos. 
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I. INTRODUCTION 



The evidence for non-zero neutrino masses and cosmological matter-antimatter asymmetry Q provides a strong 
indication for physics beyond the Standard Model (SM). An attractive way to explain these phenomena involves the 
inclusion of heavy right-handed (RH) Majorana neutrinos to the leptonic sector. As a result, tiny neutrino massess 
may be generated via the see-saw mechanism Q while the cosmic baryon asymmetry can be explained by thermal 
leptogenesis @. 



One particularly beautiful aspect of leptogenesis is that it offers a link between the baryon asymmetry and SM 
f***. i neutrino properties. The mechanism works by assuming that heavy RH neutrinos interact with the ordinary light 
neutrinos in a CP and lepton number (L) violating way in the early universe, and hence, generating an L asymmetry 
which, given the right conditions (100 < T < 10 12 GeV) @, would eventually be partially converted to a baryon 
O 1 number (B) asymmetry by non-perturbative sphaleron processes. In addition, through the see-saw formula, relation 
of the mass scales of the light and heavy neutrino spectra can be established. Therefore, the requirement of successful 
J^" 1, leptogenesis naturally leads to limits on neutrino masses and mixings. 

43 : 

Given this intricate connection, it is of great interest to explore possible neutrino models that can give rise to 
successful leptogenesis. By the same token, it is intriguing to ask what are the implications of a given model, which is 
specifically designed to address a different issue, in the leptogenesis context. The aim of this paper is to conduct such 
' studies on a specific class of see-saw neutrino models (as proposed by Low [l[) that has an Abelian family symmetry 
and an extra real singlet in the Higgs sector. Because of the interplay between the new ingredients, these models 
predict a fully hierarchical light neutrino spectrum (ie. mi = 0), as well as, #13 — in the mixing matrix. More 
importantly, they contain fewer free parameters than the default see-saw, which is the original motivation behind the 
previous work of Low. The analysis here will focus on the suitability of these models for leptogenesis in several of the 
typical setups. 



Over the years, there has been a dramatic increase in the sophistication of the quantitative analysis of leptogenesis. 
Many previously neglected effects such as thermal corrections B, different washout [1, [j| E3, EH EH an d spectator 
processes [H, E3| and, above all, flavor effects [HI, EE Ell EH [190201 12H H3| have been considered in recent analyses. 
Other variations to the general scheme, including resonant leptogenesis [23| . asymmetry production dominated 
by the decays of the second lightest RH neutrinos [2J| and models with more than three heavy RH neutrinos 
[25(, have also received attention. However, for definiteness, in this paper we concentrate on the case of two or 
three heavy RH neutrinos (depending on which sub-class of Low's models we are considering) and that their mass 
spectrum is hierarchical, ie. Mi <C M2 <C M3. Furthermore, we assume that the final L, and hence, B asymmetry 
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is predominantly due to the decays of Ni's. The present investigation makes use of various existing results and 
attempts to draw comparsions where appropriate with the hope to deduce whether these specific neutrino models 
may favor certain regimes or predict significant deviation from the standard scenarios (for both the one-flavor and 
the multi- flavor cases). In most situations, we work in the "strong washout" regime although some discussions on 
the "weak washout" regime are also included. 

With these motivations in mind, the paper is structured as follows. Section HT1 provides an outline of the neutrino 
models under consideration. It serves to highlight all the key features, as well as to recast them in a form suitable 
for later discussions. We also comment on the fine-tuning of parameters in these models. Section UTT1 begins with a 
brief overview of the leptogenesis scenario in the one-flavor case. The aim is to highlight the important issues and 
parameters so that a comparsion between the default sec-saw and our specific models can be made. In Section IIV1 
the analysis is repeated for the case when flavor effects are included. A summary of results is presented in Section [V] 



II. SEE-SAW MODELS WITH ABELIAN FAMILY SYMMETRY AND AN EXTRA REAL HIGGS 

SINGLET 

While the SM has been very successful in explaining the dynamics of sub-atomic particles, it is not without its 
shortcomings. A typical example is its inability to predict quark and lepton masses. Such parameters are put in by 
hand using data from experiments. Altogether, there are almost 20 free parameters in the minimal SM; and should 
neutrinos gain mass via the see-saw mechanism, even more parameters will be required. Therefore, from the model 
building point of view, it is natural to look for ways to reduce the number of parameters. 

A first step towards this may be achieved by enlarging the symmetry in the leptonic sector so that the symmetry 
group becomes 

G = SU(3) C <g> SU(2) L <g> U(1) Y x G 

family j (i) 

where Gf am ii y is a leptonic family symmetry. But it has been shown in [26j that models with an unbroken family 
symmetry and just the SM Higgs do not have more predictive powers than the default see-saw. Hence, one would 
need to expand the Higgs sector at the same time. 1 The exact nature of the coupling between the lepton and Higgs 
fields dictates the neutrino mixing matrix, C/pmnSi produced via family symmetries. In the work of Low [U.[2^|. effects 
of Abelian family symmetries with additional Higgs singlets, doublets or triplets on certain features of the mixing 
matrix were studied. It was found that the simplest models that can predict #13 = contain one extra real Higgs 
singlet which transforms non-trivially in family space. 2 The action of the family symmetry with one singlet can be 
summarised as follows. Suppose the leptonic Yukawa and mass terms are: 3 



£ma SS = —L Y\ <f> l R -LY v <S>v R - - {v R Y Y xX vr~ ~{vr) c M baro v R + h.c, (2) 

where Yi,Y u ,Y x and Mb are are coupling matrices in the family basis. L = (vl,Il) T is the left-handed (LH) lepton 
doublet while Ir and v R are the RH charged and neutral lepton singlets respectively. The Higgs fields are $ (SM 
complex doublet) and \ (real singlet) with $ = ia 2 <$>* . Family symmetry, Gf am ii y , demands that the full Lagrangian 
is invariant under the unitary transformations, 

L^SlL, Ir-^Si r Ir, vr^S Ur vr, $ — * 5$$ and x ~ * S x x, (3) 

in family space. These transformations restrict the coupling matrices to certain forms and thus it is possible to 
generate a mixing matrix with U e3 = 0. Here, J/pmns is parametrised by [3(| 

(C12C13 S12C13 si 3 e~ i,5 \ 

-S12C23 - C12S23S13 e lS c 12 c 23 - s 12 s 23 s 13 e lS s 23 c 13 x D v , (4) 
S12S23 - C12C23S13 e lS -c 12 s 23 - S12C23S13 e lS C23C13 / 



1 Alternatively, other symmetry breaking mechanisms could be used instead. 

2 The CHOOZ experiment 28] indicates that 613 is very small and, at best fit, it is taken to be zero. 

3 Note that it may be tempting to exclude the Mbare term in (f2j) and have RH neutrino massess generated by the \ term only. But it 
turns out that, without Mt, are , the solar mixing angle, 612, is forced to be maximal by the symmetry (which is incompatible with the 
best fit data [2!). 
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TABLE I: Diagonal representation of the transformations in family space that gives U e 3 = for three or two (in Case 3) RH 
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TABLE II: Coupling matrices generated by the transformations in Table U where "x" denotes an arbitrary complex entry. Note 
that Y x and Mbaro are symmetric matrices. 



where D v = diag(e' ¥ ' 1 / 2 , e l(p2 ^ 2 , 1) is the matrix containing all the Majorana phases while S is the CP violating Dirac 
phase, and s mn — sm.6 mn ,c mn — cos0 mn - Therefore, U e 3 — corresponds to 6*13 = in this parametrisation. A 
possible representation of the sets of transformations that leads to this is shown in Table |TJ Their associated Abclian 
symmetry and coupling matrices are presented in Table [TTJ Although the analysis is done assuming that there are 
three RH neutrinos, the desired property of C/ e 3 = can also be generated by models with only two RH neutrinos. 
This is evident from the vanishing third column of Y v in Case 1 of Table HU Hence, in this case, the third RH 
neutrino is actually decoupled from the LH sector. If we remove it from the theory, Case 1 is reduced to Case 3 (see 
Table III and ITT|) , 4 It should be noted that these models keep the atmospheric mixing angle, #23, and the solar mixing 
angle, 12 , as free parameters. Despite this and the more complicated Higgs sector, it has been shown that, overall, 
these models have less arbitrary parameters than the standard see-saw. 



Furthermore, the texture zeros in the coupling matrices together with the see-saw formula 



m v ~ -mjjMjj mo, (5) 



give rise to an m v (in any basis choice for mi and Mr) which has the following form: 

/ ai a 2 bi a 2 &2 

a 2 6i a 3 bj a 3 bib 2 I , where a±, a 2 , a 3 , 61, b 2 G C. (6) 



ya 2 6 2 o 3 6ife 2 a 3 b\ 



In the above, fhn — Y v {<&),mi — and Mr = Y x {x) + Mbare where ($) and (x) are the VEV of fields <& and x 

respectively. This structure of m v in (J6j) predicts one of the Majorana mass eigenstates iv) of the light neutrinos to 



4 In fact the physics of Case 1 and 3 are equivalent, and henceforth, we shall not discuss Case 3. 
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TABLE III: The structure of the coupling matrices: Y v in the family basis and h v in the mass eigenbasis of the charged leptons 
and RH Majorana neutrinos. In general, "x" denotes an arbitrary complex number. 



be massless. For better illustration and subsequent discussions, it is convenient to rewrite the Lagrangian of {2J in 
the mass eigenbasis of the charged leptons and heavy RH Majorana neutrinos: 

Anass = --ND M N-~lh t $e-Th v $ N + h.c, (7) 

where I = (v' L , l' L ) T and e are the charged lepton doublet and singlet respectively. We have defined the heavy Majorana 
neutrino field: 5 N = {v' R + v'£)/y/2. The charged lepton mass matrix is given by D\ = diag(m e , m M , m T ) = hi{<&) 
while for the heavy neutrinos, the mass matrix is Dm — diag(Mi, M2, M3). Although m v produced via the see-saw 
formula ([5]) is unaffected by the basis change, in general, h v in J7]) would have texture zeros different from Y v (in fact, 
the texture zeros would disappear in most cases; see Table HIT)) . Using U from (|3]) to diagonalise m v and choosing the 
sign convention that D m = —U^m v U*, one gets D m = diag(mi, rri2, m.3) with mi = and mi < TO2 < 1TI3. Therefore, 
the light neutrino mass spectrum is fully hierarchical in these models. Relations between the m^'s can be established 
by invoking the neutrino mixing data: 6 

m\ - m 2 2 = Am 2 tm (Am 2 sol ), (8) 

m 2 ~ m l = Am L ( Am atm - Am sol)- ( 9 ) 

When mi = 0, one then gets 

m l = m ltm ( Am ltm), (10) 

ml = Am s 2 ol (Aml tm - Am s 2 ol ), (11) 

where m atm = J Am 2 tm + Am 2 ol « 0.05 eV. 7 It should be noted that the conclusion in this section remains valid 
under one-loop renormalisation group running [l[. 



A. Parameters fine-tuning 



In order to make these models workable, the Higgs sector was expanded to accommodate the real singlet %• H s 
addition has inevitably introduced a new energy scale, (x), to the theory. Since Mr = Y x (x) + A/bare depends on 
this, it is essential to understand the implications of the scale of (%) in relation to other parameters in the model. 
The most general and renormalisable Higgs potential incorporating x is given by 

V{*, X) = \ti*{&*) + jA ($t$) 2 + I /ixX 2 + 1a xX 4 + Hx^^)x\ (12) 

where jij, and fi x are in general functions of temperature (T). Also, if the potential is to be bounded from below, 
then /^0 X > — ^J\^\ x j2. From the see-saw mechanism, we expect that (x) S> ($), and hence, T CiX ^> T c ^ ~ O(10 2 ) 
GeV, where T c denotes the critical temperature for symmetry restoration. In fact, one would want T c . x > T xe h, the 



5 In general, v' L ^ v (the mass eigenstate for light neutrinos). v' L and v' R are new fields from the change of basis. 

6 The inverted hierarchy case is in brackets. 

7 Aml tm = (2.6 ± 0.4) x 10" 3 eV 2 [UH and Am s 2 ol = (8.0±°;|) x 10" 5 eV 2 HHH- 
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reheating temperature, so that any topological defects (domain walls) created by the spontaneous breaking of the Z2 
discrete symmetry of x are eliminated via inflation. The required hierarchy, (x) 3> ($), is ensured (to tree-level) if 
M^x — * 0- This can be seen from the tree- level minimum condition for (fT2|) : 

M0 = -A ($) 2 - 2^ X { X )\ (13) 
^ = -A x ( X ) 2 -2/i 0x ($) 2 . (14) 

If 0(X c f > ) ~ C(A X ) = 0(1), then /i^ x — > guarantees that remains at 0((&) 2 ). In this limit, it also means that <I> 
and x fields are decoupled from each other. 

Since, in a typical thermal leptogenesis analysis it is usually assumed that the reheating temperature (T rc h) after 
inflation is larger than the mass of the decaying heavy neutrino, 8 altogether, one has: 

T c , x > T rch > Mi or 2 » T C)0 . (15) 

This relation implies that Y x must be suitably fine-tuned so that Y x (x) — O(Mj). In other words, elements of Y x 
cannot be arbitrarily small, and as a result, radiative corrections to the potential V($, x) due to x's coupling to other 
fields in the model may destroy the hierarchy between ($) and (x). Nonetheless, for the purpose of this paper, we 
shall work in the assumption that such stability problem and any additional fine-tuning of the framework will be 
inconsequential to our main discussion. 



III. IMPLICATIONS IN LEPTOGENESIS WITHIN THE ONE-FLAVOR APPROXIMATION 



In order to understand the potential implications of our specific models on leptogenesis predictions, it is essential 
to get oriented in the standard scenario. The key relation that captures the dependence of the predicted baryon to 
photon number ratio (j]b) at recombination time on the elements of thermal leptogenesis can be written as [34j 

3 

~ de\ k[ , (iVi-dominated case) (16) 

where d is the dilution factor that accounts for the partial conversion of the generated i]b-l into t\b through sphaleron 
processes, as well as the increase of photon number per comoving volume from the onset of leptogenesis to recombi- 
nation; Ej measures the CP asymmetry in the decays of Nj: 

e _ TjN, -> 1$) - r(jv,- -> m _ T 3 - T 3 

3 T(Nj Z$) + IX^- — >• I$t) Tj+Ti ' 

while represents the (final) efficiency factor for B — L production from Nj decays which takes into account the 
initial conditions and the dynamics of particle interactions in the leptogenesis era, in particular, the interplay between 
decays, inverse decays and AL ^ scatterings. Hence, there are potentially three places where our models may 
modify the overall t\b prediction. 



A. Dilution factor 



To track the time evolution of the number density of a quantity, X (eg. Nj, B — L01B), in an expanding universe, it 
is more convenient to consider the number of particles, J\fx(t), in a portion of comoving volume (R 3 (t)) that contains 
one photon at time t' <§C t\ op to 9 than the conventional, nx(t). We choose the normalization for R 3 (t) such that in 



This is done so to minimise theoretical uncertainties although in certain regimes, this restriction may be relaxed without any appreciable 
change to the predictions. Here, the decaying neutrino refers to the one that dominates the L asymmetry generation (usually Ni). Note 
that in the hierarchical limit for RH neutrinos, N3 decays cannot be the dominant contribution to the final asymmetry |33|l . 
ilepto denotes the time at the onset of leptogenesis. 
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relativistic thermal equilibrium, it contains on average M^. (t <C ii ep to) 
asymmetry at recombination time, we have 



Vb = 



M% 



■A/^y (tree ) 



1 heavy RH neutrino. Hence, for the baryon 



(18) 



where M B denotes the final B value after leptogenesis and sphalerons. In the simplest case of constant entropy (s) 
and assuming standard photons production from t' to i rcc , one has 



A/" 7 (t roc ) 



M^(t') 



4 x g* a (t') 

3 g* s (tree) 



4 434/4 
3 X 43/11 



37. 



(19) 



The pre- factor of 4/3 originates from our choice of normalisation. g*(t) is the relativistic degrees of freedom at 
time t. In (|19|) . we have already assumed the ^-dominated scenario with g*(t') taken to be the degrees of freedom 
from all SM particles and Ni only. Up to this point, the calculation is as per usual because the only new parti- 
cle in our models is the physical H x which gains mass at a very high energy (T CiX ), and at time t' , it is non-relativistic 
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Another element that contributes to the dilution factor comes from the imperfect conversion of M B _ L intoA/'Jj. The 
simplified picture which assumes that sphaleron processes are only active after the leptogenesis era, gives a sphaleron 
conversion factor of a spn = M B /M B _ L = 28/79 (3f| in the SM with one Higgs doublet. 11 In a more thorough analysis 
when the combined effect of all spectator processes [ill, [l4| in the plasma (eg. Yukawa interactions, QCD and 
electroweak sphalerons) is taken into account, the resultant value receives a 20% to 40% enchancement or suppression 
[l4j depending on the specific leptogenesis temperature Ti op to assumed. 12 This change comes about because chemical 
potentials of particles in thermal equilibirum are modified during leptogenesis, and not all of these potentials are 
independent due to SM, Yukawa and sphaleron interactions [35|. As a result, washout processes 13 that control 
M B _ L get enhanced or suppressed, leading to an overall change in a sp h- The dependence on temperature is origi- 
nated from the fact that an increasing number of Yukawa or sphaleron processes comes into equilibrium as T decreases. 



Applying these ideas to our models, it is not hard to see that a sp h receives no modification from the existence of H x . 
This is because none of these interactions: H x «-> NjNj, NjNj <-* NkNk (via H x ) changes Ml or M$ in the plasma, 
and for most temperatures, they are not in equilibrium anyway due to the heaviness of H x and Nj. Furthermore, even 
if H x is light, the process H X H X <-> $ $t which can change A^ (and hence Ml) is impotent since /j^ x — > 0. Therefore, 
our specific models do not change the prediction of the dilution factor d with respect to the standard see-saw setup 
no matter whether spectator processes are included. 



B. CP asymmetry 

While dilution d and efficiency Kj (see Sec. lIII Ct govern the portion of the generated asymmetry that would survive 
after the entire process, it is the CP and L violating decays of the heavy neutrino (Nj — » or 7$t) that give rise 
to such asymmetry in the first place. This quantity is defined as in (|17[) and its value is controlled by the so-called 
see-saw geometry [33j. Note that at tree level 

r . =r . = M^k M . ; (20) 



In general, Mjj x depends on A x which is unknown. But as in Sec. Ill Al we assume that 0(\ x ) ~ O(A^) = O(l), which is not 
unreasonable if one expects to find the Higgs (H^) at TeV scale. As a result, this implies that Mjj x ~ C((x)) i s very heavy. Moreover, 
even if H x is included in the calculation of ffj(t'), it will only change J 19D by less than 2%. 

1 This factor is somewhat different if electroweak sphalerons remain in equilibrium until slightly after T c $ [36| . Note that this case is 
highly probable because the electroweak phase transition seems to be not strongly first order. 

2 The temperature region which Ti epto falls into also dictates the importance of flavors in leptogenesis. We defer our comments on this to 
SecQV] 

3 These are the same washout processes that determine the efficiency factor . We include these effects here rather than in Sec. IIII CI 
because they are related to spharlerons. But in cases when flavor effects are important, it would be more natural to incorporate them 
in the discussion of «J a (the flavor dependent efficiency factor). 
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where we have used the notation of ([7]). Thus, no asymmetry is generated at tree level. The leading contribution 
to Sj comes from the interference of the one-loop vertex and self-energy corrections with the tree level coupling. A 
perturbative calculation of this yields [33, [ijH 



Im 



{hlK)l 



(hlh„) 



Mi 

Ml 



+ fs 



Ml 
M l 



where fv{x) and fs{x) are given by 

fv(x) = ^ 



1 - (1 + x) In 



1 + x 



and f s (x) 



1-x 



(21) 



(22) 



which denote the vertex and self-energy contributions respectively. 

For j — 1 and in the limit of hierarchical RH neutrino with M\ <C M^^, (ie. i > 1), we have 



fv(x)+f s (x) 



2Vi 



Therefore, the CP asymmetry for N\ decays in this limit is given by 



3Mi 



l&ir{hlh v )ii h , x 

3Afi 
l&ir{hlh v )ii 



Im [(hlKD^hlK)^] . 



(23) 

(24) 
(25) 



It is customary to rewrite (|25p in terms of a complex orthogonal matrix (Q) whose relation to the Yukawa coupling 
matrix h v is given by [HJ 



n = (*) d-^i 2 k Dm /2 = D7 n 1/2 m D d~; /2 with n T n = i, 



(26) 



where all symbols are defined as in Sec. [II] with matrix B = A^ 1 / 2 defined by B 2 = A^ 1 and AA~ X = /. In deriving 
([2"6"| . we have used the see-saw relation (JS|) in the basis where the charged lepton and heavy neutrino mass matrices 
are real and diagonal (ie. m v ~ — mD-D^m^) and the definition of D m = —U^m v U*. Using ([2"6"| in (|25p and after 
some manipulations, one gets [39[ 



3Mi E fc ^M^i) 

3Afim atlr 
16tt($) 2 



f3(mi,mi,n 2 kl ) , 



(27) 
(28) 



where we have introduced: the effective neutrino mass, rri\ and a dimensionless quantity, /3(mi, mi, QjL), which are 
defined as [911331] 



mi 



(m^TO D )n 



and 



/3(mi,mi,fifci) = 



Mi 



E fc m 2 fc Im(Q 2 fcl ) 
m a tm Efc mfc If^J 
E fc ^Im(fi 2 fel ) 



m a tm mi 

respectively. Note that one can easily get to ([31]) by using relation (|26| in ([29 



(29) 

(30) 
(31) 



It will become apparent later (see Sec. MI C[) that the value of mi plays a crucial role in governing the washout 
regime that is relevant, as well as determining the final efficiency factor. Suppose that rh\ and M% have been fixed due 
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to other considerations (eg. washout regime selected or potency of flavor effects), expression ([28]) then implies that 
the size of the CP asymmetry is controlled only by mj and the configuration of the three S^'s (which are all neutrino 
model dependent parameters). Although these f2| x 's are generally complex, the orthogonality condition (^ fe Vl\ x = 1) 
means that only 3 real independent parameters are needed to specify them. Moreover, since the numerator of (|31[) 
depends on the imaginary part of 0| 1; only 2 of these 3 parameters will manifest itself in /3(mi, mi, ^l^i) wnen 
mi = Sfc m feMfeil i s a constant. This can be made clear if we define f^i = Xk + iYk and use the orthogonality 
condition to re-write (|31[) as 

o, ~ n2 \ Y 2 {m 2 - m 2 ) + Y 3 {m 2 -m 2 ) 

0{m 1 ,m 1 ,tt kl ) = . (32) 

wi atm mi 

An important question to ask though is which sets of fij^'s will yield a maximum CP asymmetry given fixed values 
for mi and m\. By the same token, it is useful to ascertain the upper bound 14 for 0{m\, mi, f2|i) for a specific class 
of neutrino models. For this purpose, it is convenient to break up the dependence on mi, m,i and Q|i and introduce 
an effective leptogenesis phase Sl (3a, liol. [ill |42| 

P(rni,mi^ 2 kl ) = /J max (mi, mi) sin ^(mi, m-i, fi^i), (33) 
= 0i(m 1 )0 2 {m 1 ,m 1 )shi6 L (m 1 ,m 1 ,Q, 2 kl ), (34) 

where /3 max < 1 represents the maximal value for given a particular mi and rh\. Clearly, the CP asymmetry is 
largest when the configuration of fi^i leads to sin<5i = 1. The upper bound for when only mi is fixed is denoted 
by function 0i while 02 is a correction to 0i when mi is also fixed at a given finite value. By analysing (f3"2"|) , it can 
be shown that [H, S3 

/MmiH^^i. (35) 
m atm 

To get 02, one observes that for a generic f2 matrix, a configuration that maximises while keeping mi the same is 
achieved when Sl^ = [33j. This allows one to rewrite (|3^|) as 

y,(mg-mj) 
m a tm mi 
_ Y 3m sm5 L {ml - ml) 
m atm mi 

where Y 3m is the maximum value of Y 3 when Q21 = 0- By putting (|3"4"|) . ([55]) and (|3T|) together, one gets 

0i = --^- = ^^Y 3m . (38) 
01 smfli mi 

In general, Y 3m will depend on the light neutrino masses. For a fully hierarchical neutrino spectrum (mi = 0) like 
our models with family symmetry, we have 0\ = m 3 /m atm and 2 — 1 (with X 2 = Y2 = X 3 = and Y 3m — mi/m 3 ). 

Using the above decomposition of the CP asymmetry, it is then very easy to deduce the specific predictions of 
our models using the structure of h v . Starting from the h„'s given in Table [mi 16 we first note that for Case 2, the 
corresponding f2 matrix is arbitrary since all entries of h v in this case are unconstrained by the symmetry. As a 
result, from the above general analysis, it is clear that Case 2 certainly has enough parameter freedom to create a CP 
asymmetry for the leptogenesis purpose, and in fact it predicts no additional restrictions other than those implied by 
the model independent analysis. 

The situation for Case 1 is slightly different though due to the presence of texture zeros. For the purpose of working 
out f2, we replace all texture zeros in h v by an infinitesimal parameter g. Hence, we have for this case 

(39) 



(37) 





fx 


X 


Q 


h v = 


X 


X 


Q 






X 


Q 



Without loss of generality, we may adopt the convention that the 
< PimufhunlJ. 

For normal heirarchy, 7713 = m a tm by definition and so fi\ = 1. 
We shall ignore Case 3 for it is effectively identical to Case 1. 



CP asymmetry is positive and hence the lower bound is given by 
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Applying this in (126)) yields an orthogonal matrix of the following form: 



1 

y/m 2 M 1 
1 



OC 

OC 



VmiM 2 
1 

1 



OC 
OC 



g 

M3 



(40) 



Voc — _ (x — — 1 

\/m3 Mi \fm^M^ J 



where we have kept all m^'s and Mj's as variables while assuming J7 e3 = in deriving this. In the limit g — > 0, one 
can immediately conclude that f2 2 3i^33 ~~ > because 7712,3 and -^3 are finite and non-zero. On the other hand, this 
model predicts mi = and therefore, at first glance, O13 is indeterminate while On and f2i2 are infinite as mi, g — ► 0. 
But by the orthogonality condition, both issues can be resolved and f2 simplifies to 







fi= y/T 



n 2 
"31 







(41) 



n 3 i vi 



It is important to realise that this orthogonal see-saw matrix is identical to the special form derived from models with 
A/3 — > 00 (or effectively models with only two RH neutrinos) 0, [44|. 17 The only distinction between these generic 
models and ours is that (|41[) is originated from the Yukawa couplings and family symmetry (which gives mi — as a 
by-product). Therefore, the constraints from our Case 1 on the CP asymmetry and their subsequent implications in 
thermal leptogenesis would be very similar to those models with only two RH neutrinos [45| . 



The first point we notice is that, in Case 1, mi is bounded from below. This can be illustrated by using (I41[) in the 
definition of rh\ and one obtains 



mi = m 2 |l — O31I + 7713IQ31I 



(42) 



When \Vl\ x \ < 1, then mi ~ m-2 whereas when |f2§i| S> 1, then mi 3> m 2 . Thus, we have mi > m-2. Notice that in 
the previous case, we can safely set X 2 = Y% = X3 = in order to maximise fti which leads to a relation rh\ = m^Y^, 
and hence, no bounds on rh\. But in the current case, X2 and I2 are no longer arbitrary as the freedom to change 
them while keeping fh\ constant has been lost due to the fact that flu = 0. Because of this reduction of free variables, 
it is more straight forward to consider the function (3 directly which is now modified to 



= 



(m 



m 



m atm mi 



(43) 



with 



mi = m-2 



m 3 



(44) 



where we have imposed X2 = 1 — X3 and Y2 = —Y3. As a result of the bound: mi > ma, the size of (3 (and hence the 
CP asymmetry) is very sensitive to the size of Y3. This is because, unlike before where I3 3> 1 and I3 -C 1 lead to 
rh\ 3> 1 and rh\ <C 1 respectively, which then means that (3 (oc Y^/fhi) approaches the same value (~ 1) for the two 
limiting cases, in ([4*3")) we have instead 



(m 



x < 



m 2 



Yo 



. (m 3 + m 2 )Y 3 





m 3 - m 2 
m atm 



ifr 3 «A 3 «l, 
if y 3 > A 3 > 1. 



(45) 



Therefore in order to obtain maximum CP asymmetry, mi must be very large. However, this is potentially detrimen- 
tal to the success of thermal leptogenesis because rh\ controls the potency of the washout rates (hence the size of the 
final efficiency factor, see Sec. IIII CI) and is therefore upper bounded if the correct L asymmetry is to be generated. 
Fortunately, Mi also dictates the final efficiency factor and the net result is that in order to circumvent the problem, 
one requires a larger lower bound for M\ at the same time (Tlj . Within the non-supersymmetric context, the mass 



17 This observation also justifies our claim that Case 1 and 3 are in fact equivalent. 
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of Mi, which is related to T ro h, is not tightly constrained and thus Case 1 of our model will be workable in leptogenesis. 

Another observation from (|45p is that the maximum attainable (3 value in the limit of Y3 3> A3 3> 1 is drastically 
different depending on the hierarchy scheme assumed for the light neutrinos. For normal hierarchy, /3 max « 0.82 
whereas in the inverted case /3 max ~ 0.18. Although this result alone is not sufficient to rule out the inverted case, it 
is clear that, for Case 1 at least, the inverted hierarchy is strongly disfavored. 

Overall, we have demonstrated in this subsection that our specific neutrino models with family symmetry can 
naturally generate the required raw CP asymmetry necessary for successful leptogenesis. The parameter space in 
these models is not overly restrictive, and the predictions are identical to a couple of special cases in the default 
see-saw, namely, the fully hierarchical light neutrinos limit (Case 2) and the two RH neutrinos scenario (Case 1 or 3). 



C. Efficiency factor 

The dynamical generation of a B — L asymmetry in the leptogenesis era depends on the out-of-equilibrium decays 
of the heavy Nj's, as well as other interactions in the thermal plasma. These non-equilibrium processes which control 
the evolution of fifm and Mb-l are quantified by a system of (at least) two Boltzmann kinetic equations. An 
important issue in setting up these is the identification of all relevant interactions which can modify A/jv,- and Nb-l 
(see Table HV| . In the minimal setup where decays of the lightest RH neutrino N\ dominate the final B~L asymmetry 
and assuming hierarchical RH neutrino masses, only the evolution of AV, and Mb-l are relevant and the kinetic 
equations (in the default see-saw case) have the follow form [8j, [9J, [ll|, uM, |46[ : 



dMi 



dz 
dN B -L 
dz 



-{D + S)W Nl -M£), (46) 
-ei D (Af Nl - A/^) - W Mb-l, (47) 



where the dimensionless variable z = M\/T is defined for convenience. A/^ denotes the equilibrium value for A^Vi 
which is now a function of z. The terms D = Td/(zH), S = Tg/(zH) and W = Tw/{zH) encapsulate the reaction 
rate of the various processes with H being the Hubble expansion rate which is given by 

H ~ 1.66-v/gJ — ^ — x , 48 
v 5 Mpi z 2 v ' 

where g*(z ~ 1) = 106.75 is the number of relativistic degrees of freedom, 18 Mpi w 1.22 x 10 19 GeV is the Planck 
mass. The thermally averaged (total) decay rate To which accounts for decays and inverse decays (N\ <-> Z$), is 

related to the zero-temperature decay rate T y D = T\ + T\ via 



K,2{Z) ' 



T^'^fl, (49) 



where JC n (z) is the nth order modified Bessel function of the second kind. T$ represents the AL = ±1 scatterings 
(eg. N\l «-> tq). Typically, scatterings which involve gauge bosons are ignored to first approximation. 19 The 
washout rate Tw, which incorporates everything that tends to erase the B — L asymmetry, is depended on the rates 
for inverse decay (Z$ — ► Aq), all AL = ±1 scatterings (except those with V^), as well as, the AL — ±2 processes 
mediated by Ni (eg. «-> 1&). 20 

To describe the behavior of the solutions to (|4^|) and fTT)) . it is customary to introduce the decay parameter [48[ 

p(T-0) ~ 



This includes all SM particles only. The N\ degrees of freedom is not included because in the preferred strong washout regime, they are 
non-relativistic at ti op t (corresponding to 2i epto = Mi /Ti epto ~ 1). 

It turns out that they are only critical if one considers the weak washout regime. See for example 0, l23t |47| 
AL = ±2 processes mediated by 7V2,3 are suppressed at z ~ 1 because M2, 3 Mi. 
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AL — 




AL = ±1 




NjNj (tree + loop)* 


Nj «-> 


(tree + loop) 


NjNj 


<-> N k N k (H x ,s)* 


Njl<r+tq s) 








-►!«(*,*) 


NjNj 


^NjNj (H x ,s)* 








^NjNj (H x ,t)* 






NjN k 


^NjN k {H x ,t)* 




+ * v M (*,«) 






Njl++$ V„ {I, t) 




<-> H x &$ {H x ,s)~k 




<-> $ i ($,t) 


N 3 H X 


~Nj&<f> (H x ,t)* 




<-> $ z (ii, t) 








<-> v M i t) 


AL = ±2 




" V M 1 (1, a) 


1$ * 


-> Z$ f (Nj,s) 






1$ * 


^I$t {Nj,t) 


Z$ ^ 


N 3 H X (N j>S )* 


11^ 


-> $ t $ t (N^t) 


Njl*- 








H X U 





TABLE IV: A collection of potentially important processes in leptogenesis. The type of process (eg. tree-level, vertex/self- 
energy loop, s- or i-channel) and, where applicable, the mediating particle are in brackets. Interactions that are not in the 
default see-saw are marked by a "*". q, t and denote the (3rd generation) quark doublet, top quark singlet and gauge boson 
respectively. 



where 

16 . F^tf <$> 



3 V 5 Mpi 



O(10" 3 ) eV, (51) 



is the equilibrium neutrino mass, mi is the effective neutrino mass as defined in (|29[) which measures how strongly 
coupled Ni is to the thermal plasma. It should be noted that the size of rh\ relative to m* marks the boundary 
between the so-called weak (mi < m*) and strong (rh\ > m*) washout regimes where the corresponding analyses 
are qualitatively different [49|]. This dependence of the solutions comes about because the interaction terms D, S and 
W\ (defined as W\ = W — 6W, where SW represents the contribution from non-resonant AL = ±2 processes) are 
proportional to fh\ [TlT |: 

Mpi mi „ Mpi Mi m 2 , , 

D,S,Wkx-^±, while SW cx P 1 I , (52) 

where m 2 = mf + m 2 , + m§. For hierarchical light neutrinos and Mi <C 10 14 GeV (note that both of these are 
satisfied in our models 21 ), it turns out that contribution from SW can be safely neglected [HI, EH, and therefore 
the generated B — L asymmetry is to a good approximation independent of Mi. In the interesting case of strong 
washout, the W\ term is dominated by inverse decays @, llsl- Moreover, in this regime, the heavy Nis in the plasma 
can reach thermal abundance before fi ep to even if the scattering term S is turned off. As a result, the details of iVi's 
production prior to their decays become irrelevant and all subsequent analyses are greatly simplified. In the light of 
this, we shall concentrate on the strong washout regime in much of our discussions here. 



dAf Nl 
dz 



With the above simplifications, (|46| and (1471) reduce to 

-DW^-Affi), (53) 
-ei D (Af Nl WFNb-l , (54) 



dz 



The only restriction on M\ is coming from II15I I which does not limit the possibility of Mi -C 10 14 GeV. 
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where W]°{— W\) is the dominant piece in the washout term which originates from inverse decays. 22 The solution 
for Nb-l can be expressed in an integral form [Zij, l5l| 



M B - L {z) =Af B _ L e-^ dz ' wl ° {z,) -e lKl (z), (55) 

where N"b-l an d z i denote the initial values for Mb-l and z respectively. The corresponding efficiency factor is given 
by 

K 1 (z) = - f d^e-e^'O^"), (56) 
J Zi dz' 

Since in the strong washout regime one can invoke the approximation dAf^/dz' ~ dAf^ /dz', the final value as z — > oo 
may be readily worked out and one obtains [49| 

rFTv 1 " 6 )> for^>l, (57) 

where we have expressed k\ as a function of the decay parameter K\\ and zb (3> ziepto — 1) is the value around 
which the asymptotic expansion of the z' integral in (|56p receives a dominant contribution. The temperature Tb 
that corresponds to zb (= M\/Tb) is referred to as the baryogenesis temperature. For practical purposes, (|5T[) is well 
approximated by the simple power law (49j : 



(2± 1) x 10" 



0.01eVl 11±01 



m i 



for mi > m*. (58) 



Assuming we are working in the strong washout region, we will now argue that the standard results described above 
are directly applicable to our models with family symmetry. Of the AL — or AL — ±1 new interactions originating 
from coupling to H x , all involving an external H x in the initial state (eg. H x — > NjNj) can be disregarded because 
the reaction density for H x is almost zero at z ~ 1 while the reverse processes are kinematical forbidden due to the 
heaviness of H x . 23 Furthermore, the Nj-Nk scatterings mediated by a virtual H x do not play a role because in the 
strong washout regime, the analysis is insensitive to the initial Nj abundance in the thermal plasma. Hence, there 
is no new important contribution from the list of starred (*) reactions in Table IIVI to the standard scattering and 
washout terms, and we can conclude that our models do not predict a modification to the efficiency factor k\ given 
by the default see-saw case. 



Weak washout regime 



We end this section by including a short discussion of the phenomenologies that might be essential to the analysis 
of our extended see-saw models in the context of weak washout. 24 While the efficiency factor is typically enhanced 
in this regime, one obvious drawback is that the prediction is no longer independent of initial conditions as in the 
scenario for strong washout. As a result, the Nj-Nk scatterings via a H x may be significant since these can actively 
modify M l Nl in the plasma while inverse decays are no longer strong enough to ensure thermal abundance before the 
onset of leptogenesis. 

However, upon closer inspection and assuming a hierarchial RH neutrino mass spectrum, one would not expect 
a sizable change to A/"jy . This is because in Table iTVl there is no (AL = 0) creation or annihilation process which 
involves TVi's and some lighter particles. Note that one does not need to be concerned with AL ^ processess in this 
analysis because they are in general too weak to bring iVj's into equilibrium at high temperature (Ti opt o <C T < T re h). 
Thus only processes such as NiNi *-* NkNk (k — 2, 3) can bring iVi's into equilibrium. But as has been shown in a 
similar model in 0] where the additional Nj-Nk interactions come from couplings to massive neutral gauge bosons 
(related to GUT breaking), these scatterings have very little effect on the final asymmetry prediction if one has a 



Note that in l|53| l and l|54ll . the resonant AL = ±2 contribution must be properly accounted for to avoid the un-physical asymmetry 
generation in thermal equilibrium 0, l5dH . 

It should be noted that the 2 <-* 3 reactions such as NjNj <-» H x && are in any case impotent because of \i$ x — * 0. 
For a complete review of the analysis in the weak washout regime, see for example [7I. |49H . 
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pronounced RH neutrino mass hierarchy. 25 In addition, the Vi-dominated approximation remains valid. 

It should be pointed out that all interactions involving an external H x remain unimportant. When z <C 1, processes 
that can produce Nj (eg. H x — > NjNj) are irrelevant because they are only potent before the inflationary stage (cf. 
(fl5]l) and any excess Nj created will be diluted away. On the other hand, for z ~ 1, reaction density for H x is almost 
zero as in the strong washout regime. 

Other issues that are potentially important in the weak washout regime include the effects of scattering processes 
that couple to gauge bosons 0, [23| (see Table IIV|) , as well as thermal corrections to Yd and Y$ @, S3 . These 
elements provides an additional source of theoretical uncertainty to the overall analysis and hence making the weak 
washout scenario less attractive than the strong washout case. For the purpose of this paper though, these effects 
unnecessarily complicate the analysis and therefore will be ignored. 

Although we have not shown by explicit calculations that predictions of our models will not significantly deviate 
from (or bias towards certain parameter space within) the standard case with weak washout, it is apparent from 
the above discussion that a lot of it will be highly dependent on the assumptions made on the mass of Nj , H x and 
their Yukawa couplings. Given that all of these are free parameters in our models, any deviations with respect to the 
standard see-saw case will therefore stay within the amount of these uncertainties. Hence, we can conclude that their 
predictions are effectively the same as the default case. 



IV. IMPLICATIONS IN LEPTOGENESIS WITH FLAVOR EFFECTS 

In analogy to the discussion in Sec. IHI Al on spectator processes, flavor effects become important when the 
assumption on the temperature range in which leptogenesis happens 26 is altered pj, GJ, Q3, M, GS M, EH H- 
This arises because Yukawa interactions become more prominent at temperatures T < 10 12 GeV. Whereas in 
the previous discussion the role of these interactions were merely to change the reaction densities in the Boltz- 
mann equations which led to what we referred to as a new effective sphaleron conversion factor, 27 in the context 
of flavor issues in leptogenesis, their effects on the CP asymmetry and (flavor-dependent) washout are also considered. 

Whilst in the one-flavor approximation, the lepton states, | lfj\)'s, 2S generated by the decays of Nj's are assumed 
to have evolved coherently during the leptogenesis era, for the flavor-aware situation, the presence of charged lepton 
Yukawa interactions in equilibrium essentially introduces a source of decoherence whereby the | £(j))'s are projected 
onto one of the three flavor eigenstates, | i a )'s (a = e,/i, r) with probability l a )\ 2 ■ Because of the difference in 
size between the tauon and muon Yukawa couplings, there exists two temperature thresholds: T° q and T° q , which 
govern the range where tauon and muon interactions come into equilibrium respectively. When T > T° q ~ 10 12 GeV, 
both the r- and /^-Yukawa interactions are out-of-equilibrium and flavor effects can be ignored. For the temperature 
range T° q > T > T° q ~ 10 9 GeV 0, H3, HH , only r- Yukawas are in equilibrium and one effectively has a two-flavor 
(r and a linear combination of \x and e) problem. Finally, in the case of T^ q > T, both reactions are strong enough to 
instigate a full three-flavor system. In our discussion here, we do not explicitly distinguish between the two situations 
and simply assume the three-flavor regime. Furthermore, spectator effects in the form discussed in [l3l [Tij will be 
ignored for brevity. 

To illustrate the origin of the flavor affected CP asymmetry, 29 we begin by rewriting Lagrangian ([7]) with sub-indices 

£mass = -~iVj D Mj Nj - l a (hi) a $ e a ~ J a {h v ) aj $ Nj + h.c, (59) 



The model presented in Q provides a good guide to the phenomenologies expected in ours with the exception that our models do not 
Dossess interactions that link Nj to SM particles. Consequently, our models would predict even less modifications than in the case of 

In the TVi-dominated scenario, this corresponds to T ~ M\. The moral is that if Mi is light enough then flavor effects become important. 

One may also simply interpret it as another form of modification to the washout terms. 

The subscript (j) highlights the fact that the flavor decomposition of Jyj can be different for each j. 

A proper discussion of this should be done within the density matrix framework [l5tl20l|. But for our purpose, it is enough to follow the 
more intuitive approach as in [l9l |. 
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where j = 1,2,3 and a = e,fi,T. It is then clear that the state \ and |I(f>) in (fT7|) are in general no£ a CP 
conjugate of each other owning to the complex entries in matrix h v . If we are considering the temperature ranges 
where the relevant lepton Yukawa interactions are either fully in equilibrium or out-of-equilibrium but not in bewteen 
such that state | l^y) will quickly decohere into one of the | l a ) states available, then we can effectively think of 
Tj as an incoherent sum over a of the partial decay rates Tj a = T(Nj — > l a &)- Likewise for | Itj)), \ l a ), Tj and 
Tj a = T(Nj — > 7 a &). Note here that | l a ) and | l a ) are CP conjugate states of each other. It is then straightforward 
to see that CP violation in Nj decays can manifest itself in two places: 

A. The amount of | l a ) and |7 a ) produced are not the same because 1 1^) and \ lu)) are produced at different rates 
which corresponds to Tj =/= Tj . 

B. The amount of \ l a ) and | l a ) produced are not the same because Tj a ^ Tj a (regardless of the relation between 
VjJj). 

Obviously, the second effect is only relevant when one is considering the evolution of the individual lepton flavor 
asymmetry L a . For the one-flavor approximation, one tracks the evolution of L = ^ Q L a instead and hence only the 
first effect comes into play. However, an important corollary is that when flavor effects are included, the associated 
flavored CP asymmetry, Ej a , can be non-zero even if Tj = Tj (ie. Bj — 0). 30 To properly quantify all these, it is 
convenient to introduce the flavor projectors [l9[ 

Pj a = ^ = \{ia)\L)\ 2 , (60) 

P JQ = ^ = |(7 0) |7 Q )| 2 , (j = l,2,3;a = e,//,r) (61) 

1 3 

where Tj — J2 a Tja an( ^ Tj = VL Tj a . Therefore, by definition, J2 a Pja — J2 a Pj a = The associated a flavor 
C P asymmetry is given by (l5l [l9l ] 

T — T TP —TP 

*-J« L jot _ L 3 r ja L J r 3a / fi r,x 

Eja ~ T J+ T, T.+T, ' (62) 

P 4- P P — P 

= tja + ej + tja gga ^ (63) 

^U + ^f, (64) 

where Pj a is the tree level contribution to the projector Pj a with P® a = P°j a , while SPj a = Pja — Pja is the quantity 
that characterises the CP violating effect of type B. From (|64| . it is clear that even if Sj — and hence Pj a £j 
(representing type A effect) vanishes, Sj a does not necessarily go to zero. In addition, 

£ SP ja = (Pja ~ Pja) = 1 - 1 = 0, (65) 
a a 

demonstrating that Sj = ^2 a Sj a will not depend on SPj ai and hence it is consistent with our claim that only type 
A effect is important in the one-flavor approximation. Putting these quantities in terms of parameters in (|59[) and 
setting j = 1 for the Aq-dominated scenario, we get [lj| 

pO _ p0 _ (X)al(Mal / fifi x 

r io — r ia — j—j" : : \ m ) 



This situation corresponds to having a real orthogonal matrix f! (cf. J27j) which in turn implies that CP is an exact symmetry in the 
RH neutrino sector |21| . An intriguing consequence of this situation is that leptogenesis is directly connected to the low energy CP 
violating phases in mixing matrix U |52|1 . 
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and [37 



£ lc 



77T7 : E Im { t^)" 1 ^")°*^ 



^(^^)ife + ^(^Mfci 



(68) 



with £i given by (|24p. Using these definitions, one can then calculate 5P\ a via 



31 



To derive the network of Boltzmann equations, firstly we note that the final asymmetry M B _ L would now depend 
on the evolution of all individual flavor asymmetries A/q q where Q a = B/3 — L a and M B _ L = ^ a A/q q . Secondly, 
the washout terms would be modified in the multi-flavor case because AL ^ interactions with <i> couple to state 
| Im) and not | l a ). This effect is accounted for by P\ a ~ P® a . If we assume again that the Yukawa interactions are 
either strongly in equilibrium or out-of-cquilibrium but not in the transition region, then the kinetic equations are 
greatly simplified as flavor dynamics due to the coherences may be neglected. 32 Furthermore, in analogy to Sec. lIII Cl 
a good approximation to the problem can be obtained by considering just the effect of decays and inverse decays 33 
Overall, we get (cf. ^ and (HD) [HIH 



dMi 



dz 



dz 



-D(Af Nl -Af^), 

-e la D (Af Nl - A/£j ) - PL Wr N Qo 



Like (|5"5")l , the solution to A/q„ can be expressed in an integral form: 

Mj. (*) = M' Qa e-< Z dz ' - e la Kla (z), 



where the flavor dependent efficiency factor is given by 

Kla(z) = - 



dz' 



In the strong washout regime, one then gets (cf. |57|) [22 

2 



Ki a Z B (Ki a ) 



la) ( 



1 - e 



(69) 
(70) 

(71) 
(72) 

(73) 



where K\ a = P^ a K\. Since Pf a < 1, we have K\ a < K\ and so (|73|) implies that washout is in general reduced 
because of flavor effects. As a consequence of this, one finds that the region of parameter space corresponding to the 
strong washout regime is also reduced [22j ■ 

From the above discussion, it is clear that much of the analysis done for our models in Sec. IIIII would be applicable 
to the case when flavor effects are included. Comparing (|6"9"|) and (fTD|) with (fS"3")) and ([51)) . we observe that P® a and 
eia are the only two new ingredients which contain all the additional model dependent information originating from 
flavor effects. Hence, it suffices to investigate the implications of our specific models on these two quantities. 

To begin with, we note that both P® a and e\ a depends explicitly on the mixing matrix U, unlike in the one-flavor 
approximation where everything can be written in terms of the see-saw orthogonal matrix f2. This can be seen from 
(f6"T)) and by expanding the two terms containing h v 's in (|6"5]) using (f2"6")l : 



[K)al{K)ak(hth v )lk 



MiMk 

($)4 



E 



(74) 



31 Even without explicitly calculating &P\ a , one can see from |[68) that the term proportional to drops out upon summing over 
a, since when combined with the pre-factor outside the brackets, it becomes (hlh v ) lk {hlh v ) kl which is real. This is an indication that 
this term actually corresponds to &P\ a type effect. 

32 These coherences, in the density matrix formulism [l5j . correspond to the off-diagonal terms, p a p of a density operator, p, which has 
the properties: p aa oc and p aa = Ml- In cases where these effects cannot be ignored, one must analyse the system in full. 
See for example [20l |. 

33 The real intermediate state contribution from AL = ±2 scatterings has been properly subtracted. 
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and 

m m ^IV^mkU a iU am- (75) 

' ' n E, m 

So, this illustrates that low energy CP violating phases may, in prinicple, have direct connection to leptogenesis 
(with flavor effects). However, these effects are usually masked by the complex phases in (ie. from h u ) and would 
enter directly only when f2 is real as stated before. In the context of our models with family symmetry, since the 
prediction is t/ e 3 = 0, the Dirac phase 5 in (l} does not enter into the theory at all, and hence CP violating effects 
arising from low energy parameters would entirely come from the Majorana phases 951,2- 

In terms of whether our models would give rise to a significant departure in the predictions of leptogenesis (with 
flavors), there are some general observations we can make. 34 Firstly, by examining ([66]) and ([68]) . we can safely 
conclude that for Case 2 of our models (see Table llllj) there is essentially no modifications to the standard see-saw 
scenario. This is due to the fact that all entries in h v for this case are uncontrained by the symmetry and therefore 
it can, in principle, accommodate all specific scenarios the default see-saw allows (except mi = must be obeyed). 
Secondly, using expression (|67[) for Pi a , we see that both Case 1 and 2 do not provide any predictions or restrictions 
on the projectors P-j^'s. 35 This is because Pf Q depends only on the 21- and 31-entries 36 of both U and f2 and we 
know from before that these are unconstrained by the symmetry (for Case 1, see (|4ip ). 



To study the implications of Case 1 on e\ ai we can apply m\ = 0, J7 e 3 = and £1 in ([41]) to expressions (|74|) and 
([75]) . A key feature which results is that both l[74"]) and ([75]) vanishes when k = 3, which is in accordance with the 
previous conclusion that Case 1 corresponds to the see-saw case with only two RH neutrinos. In fact, this is the only 
definite prediction from the model. Although one gets significant simplication to ([71]) and (|75p for particular values of 
k and a, the resulting expressions are still a function of unconstrained parameters, hence one has to impose additional 
conditions in order to restrict the result. For instance, when k = 2 and a — e, (|75j) simplifies to 



{KUiKUihlKU = ^I^MiM 2 ^jj^i! _ Q 2 x | _ m3fi 2 i(1 _ Q 2 ir)) (76) 

where we have used mi = 0, U e s — and (14*1"]) . Plugging this into ([fi"5]) and defining Q,^ = X 3 + iY 3 , we simply get 
(ignoring 0((M 1 /M 2 ) 2 ) terms): 

£le -- iR 3 /*i~ ^IM 2 l3. (77) 

lD7r(<£ > ) z mi 

This example illustrates that although the texture zeros can lead to partial simplification, in the end, one has to 
specify f2§i and U e 2 (in this case) in order to make any predictions about e\ e . Besides, the final asymmetry would in 
general depends on all of the £i a 's and not just E\ e . Therefore, all we may conclude is that Case 1 of these models 
will possess the phenomenologies given by the standard see-saw with two RH neutrinos. An extensive study on this 
particular situation in the context of leptogenesis with flavor effects has been performed recently in [2l| . 



So far, we have assumed the strong washout scenario exclusively. The reason for that is in the weak washout regime, 
the calculation of the final asymmetry relies on a good knowledge of the thermal history of all RH neutrinos in the 
plasma. Since our models do not provide specific information on that front and flavor effects cannot actually alter 
the phenomena we discussed in Sec. 1111 CI therefore, our conclusion on this regard would be similar. In the light of 
all these, it is clear that implications of our models in leptogenesis with flavor effects are basically the same as those 
from the normal see-saw (with mi = 0). 



V. SUMMARY OF RESULTS 



It is fascinating that two seemingly unrelated problems — the tiny masses of light neutrino and the cosmic baryon 
asymmetry — may be explained in a consistent manner by the mere introduction of heavy RH neutrinos to the SM. 



Note that there is nothing in all cases of our models which specifically indicates whether flavor effects should be included or not for M\ 
is a free parameter in the theory. Thus, leptogenesis with flavor effects is phenomenologically not excluded by our models. 
Again Case 3 is automatically accounted for when considering Case 1. 

No dependence on the 11-entry because mi = and contribution to the sum in J6TH automatically disappears. 
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While the former may be explained by the see-saw mechanism, thermal leptogenesis provides an attractive solution 
to the later. Given this connection between the two, it is natural to ask whether a specific class of constrained 
see-saw models which has less free parameters than the default setup, can give rise to successful leptogenesis. In this 
paper, we have performed a thorough check in this context on such constrained models as those proposed by Low [l|. 
These models obey certain Abelian family symmetry and contain one additional singlet in the Higgs sector such that 
the total number of arbitrary parameters in the see-saw theory is reduced. In particular, they predict that #13 = 
and a fully hierarchical light neutrino spectrum (mi = 0). 

By dissecting the dependence of the final asymmetry on sphalerons, CP violating decays and washouts, we 
have identified the key elements that can modify the leptogenesis predictions. Consequently, the implications of 
different see-saw neutrino models can be easily compared. It was found that in all cases of our models, leptogenesis 
predictions are almost identical to those allowed by the default see-saw model, with the exception that our model 
naturally selects the hierarchical light neutrino solution (for mi = is one of its main features). This conclusion is 
true respectively for both the one-flavor approximation, as well as when flavor effects are included since all of the 
essential elements that can change the final asymmetry turn out to be unconstrained in our models. In one case, 
the phenomenologies are very much the same as the standard situation with three RH neutrinos whereas in the two 
other cases, they correspond to the scenario with only two RH neutrinos. Furthermore, for leptogenesis with flavor 
effects, we have found that Majorana phases in the light neutrino mixing matrix can play an important role since in 
our models the Dirac phase disappears due to the fact that $13 = 0. 

It should be noted that the entire investigation has been done assuming the limit of hierarchical RH neutrinos. 
Also, we have adopted the Ai-dominated scenario in most situations. Although in much of this paper the analyses 
are done in the strong washout regime, there are clear indications that our general conclusion can be extended to the 
weak washout case. A more precise calculation on this front, however, cannot be carried out since our models do not 
provide any specific predictions on the RH neutrino sector. On one hand, this means that nothing new is predicted 
by our models, but on the other, this non-existence of heavy constraints ensures that these models lead to successful 
leptogenesis in most scenarios. 

Note added:— While this work was in preparation, it was found in [53[ that quantum Zeno effects can affect the 
conditions in which flavor effects in leptogenesis are applicable. Although this is an important observation, it will not 
alter equations (|69p and (|70|) . Hence, our analysis in Sec. HVI and its conclusion remain the same. 
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